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We propose a model to describe the statistical properties of wave scattering through a classically 
chaotic cavity in the presence of surface absorption. Experimentally, surface absorption could be 
realized by attaching an "absorbing patch" to the inner wall of the cavity. In our model, the cavity 
is connected to the outside by a waveguide with N open modes (or channels) , while an experimental 
patch is simulated by an "absorbing mirror" attached to the inside wall of the cavity; the mirror, 
consisting of a waveguide that supports N a channels, with absorption inside and a perfectly reflecting 
wall at its end, is described by a subunitary scattering matrix S a - The number of channels N a , as 
a measure of the geometric cross section of the mirror, and the lack of unitarity P a = ljv a — S\S a , 
as a measure of absorption, are under our control: these parameters have an important physical 
significance for real experiments. The absorption strength in the cavity is quantified by y a = trP a . 
The statistical distribution of the resulting S matrix for N — 1 open channel and only one absorbing 
channel, N a = 1, is solved analytically for the orthogonal and unitary universality classes, /3 = 1 and 
P = 2, respectively, and the results are compared with those arising from numerical simulations. The 
relation with other models existing in the literature, in some of which absorption has a volumetric 
character, is also studied. 

PACS numbers: 73.23.-b, 03.65.Nk, 42.25.Bs,47.52.+j 



I. INTRODUCTION 

Wave scattering experiments with microwave 0,0 and 
acoustic resonators J3| represent a fruitful field for the 
verification of Random Matrix Theory (RMT) predic- 
tions d 0- El since, in that experimental domain, exter- 
nal parameters are particularly easy to control. However, 
the cost to be paid is the presence of absorption due to 
power loss in the walls of the device used in the exper- 
iments. Since the appearance of the paper by Doron et 
al |2j which showed the drastic influence of absorption, 
many investigations, both experimental and theoretical, 
have been devoted to the study of absorption effects on 
transport properties [E S E3, E3, E3, EE EE E3, 
E El |H lUf E3 of classically chaotic cavities. 

A simple model to describe the statistical properties of 
cavities including absorption was proposed by Kogan et 
al. : it describes the system through a subunitary scat- 
tering matrix S, whose statistical distribution satisfies 
a maximum information-entropy criterion. The model 
turns out to be valid only in the strong-absorption limit. 
An alternative model proposed by Lewenkopf et al. [Tol ] 
simulates absorption by means of N p equivalent "para- 
sitic channels" , not directly accessible to experiment (in 
addition to the N physical channels) , each one having an 
imperfect coupling to the cavity described by the trans- 
mission coefficient T p . These parasitic channels can be 
interpreted in terms of the voltage-probe model origi- 
nally proposed by Biittiker [23L |2J|, where a fictitious 
lead with N p channels is attached to the cavity, each one 
with a coupling T p @. The total scattering matrix S 
for this problem is unitary and has dimension N + N p ; 
the iV-dimensional submatrix S thereof is the physical 



scattering matrix which is of course subunitary. In the 
limit N p — > 00 and T p — > 0, while the product 7 P = N p T p 
is kept fixed and interpreted as the absorption strength, 
this model was shown by Brouwer and Beenakker |8( to 
describe volume absorption, in the sense that the problem 
described by S is equivalent to one in which all the en- 
ergy levels of the closed cavity acquire a fixed imaginary 
part, associated precisely with the absorption strength. 
Although the above quantity 7 is not directly under ex- 
perimental control, it can be chosen so as to fit the ex- 
perimental data [2fjll2l| . 

Alternatively, instead of considering a cavity that ex- 
hibits volume absorption, we may think of an experimen- 
tal situation in which "absorbing patches" are literally 
attached to the inner wall of the cavity in a controllable 
fashion, giving rise to what we shall call surface absorp- 
tion. The physical properties of these patches could be 
determined and controlled independently of the cavity 
they are used with. 

In the present paper we propose a simple model to 
describe the statistical properties of wave scattering 
through a ballistic chaotic cavity in the presence of sur- 
face absorption as defined above. In the model [see Fig. 

the cavity is connected to the outside by a waveguide 
with N open modes (or channels), while an experimental 
patch is simulated by attaching to the inside wall of the 
cavity an "absorbing mirror" , consisting of a "frustrated 
waveguide" that supports N a channels, with absorption 
inside and a perfectly reflecting wall at its end: the mir- 
ror is described by a subunitary scattering matrix S a . 
The number of channels N a , as a measure of the geo- 
metric cross section of the mirror, and the lack of uni- 
tarity P a — ljv a — S^Sa, as a measure of absorption, are 
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FIG. 1: A ballistic chaotic cavity described by the scatter- 
ing matrix So, connected to a waveguide with N propagating 
modes and to an "absorbing mirror". The latter is formed 
by a "frustrated waveguide" that supports N a channels and 
contains an absorbing barrier inside and a perfectly reflecting 
mirror at its end; the absorbing mirror is described by the 
subunitary scattering matrix So- 
under our control: these parameters have an important 
physical significance for real experiments. The absorp- 
tion strength in the cavity is quantified by 7 a = trP a . 
Alternatively, we may think of attaching to the inner 
wall of the cavity N a absorbing patches, each one with 
a geometric cross section that supports only one channel 
with a known absorption. 

The paper is organized as follows. In the next section 
we study the model for surface absorption that was de- 
scribed in the previous paragraph. In Sect. IIIII we find 
the relation between the present model and the parasitic- 
channel model also described above: we show the equiv- 
alence ot the two models when, in the latter, the limit 
mentioned earlier leading to volume absorption is not 
taken and the parameters are chosen properly. In Sect. 
IIVI we show that the model can be solved analytically 
in the one-channel case and for one absorbing channel. 
Finally, we present our conclusions in Sec. E] 



II. THE MODEL FOR SURFACE ABSORPTION 

The system that we are interested in and that we shall 
study in the present section was described in the Intro- 
duction and is shown in Fig. ^ 

The scattering problem can be studied in terms of an 
N x N scattering matrix S, which in the stationary case 
relates the outgoing-wave to the incoming-wave ampli- 
tudes 0, |25| . As shown schematically in Fig. ^ S can 
be seen as the combination of a scattering matrix So 
which describes a ballistic cavity connected to two non- 
absorbing waveguides (one with N channels on the left 
and the other with N a channels on the right) and the 
scattering matrix S a that describes the absorbing mir- 
ror; S a is a fixed N a x N a subunitary matrix which can 
be considered as an input to the problem. The lack of 
unitarity of S is defined by P a = %-N a — S^S a : in the 
most general situation, S a , and hence P a , is not a di- 
agonal matrix. The overall absorption strength of the 
absorbing mirror can be quantified by 7 a = trP a . 



The scattering matrix So is of dimension N + N a and 
has the structure 

s °=(<°'0' (2 - i) 

where r is the N x N reflection matrix for incidence 
on the left, r' Q is the N a x N a reflection matrix for inci- 
dence on the right and t and t' are the corresponding 
transmission matrices of dimensions N a x N and N x N a , 
respectively. 

The matrix So will be considered to belong to one of 
the basic symmetry classes introduced by Dyson in Quan- 
tum Mechanics [261 , which we briefly recall in order to 
make a more precise statement. In the "unitary" case, 
also denoted by (3 — 2, the only restriction on So is uni- 
tarity, due to the physical requirement of flux conserva- 
tion. In the "orthogonal" case (/3 = 1) So is also sym- 
metric because of either time-reversal invariance (TRI) 
and integral spin, or TRI, half-integral spin and rota- 
tional symmetry. In the "symplectic case (/3 = 4), So is 
self-dual because of TRI with half-integral spin and no 
rotational symmetry. In the scattering problem of scalar 
classical waves, the orthogonal case is the physically rel- 
evant one. However, we shall consider below both [3=1 
and /3 = 2, in the understanding that the unitary case 
is to be considered as a reference problem, as it is often 
simpler to treat mathematically than the orthogonal one. 

We assume that the barrier inside the absorbing mir- 
ror is sufficiently far apart from both the entrance to the 
cavity at one end and the perfect reflector at the other 
end that the evanescent modes do not play a role Q. 
Then, following the combination rule of scattering ma- 
trices, the total S matrix which describes the cavity with 
absorption is a subunitary matrix given by 

S = ro+t' - 1 -^- T S a t , (2.2) 

where 1„ stands for the unit matrix of dimension n. 

The physics of Eq. 1)2. 2J1 is clear. The first term on 
the right hand side is the reflected part of the wave that 
enters the cavity; the transmitted part reaches the ab- 
sorbing mirror with the (matrix) amplitude to and the 
fraction S a (again a matrix) thereof returns to the cav- 
ity; the wave then suffers a multiple scattering process 
between the cavity and the absorbing mirror, until fi- 
nally the fraction t' leaves the cavity through the lead 
on the left. 

The statistical features of the quantum mechanical 
scattering produced by the cavity shown in the lower 
portion of Fig. ^ -assumed to have a chaotic classical 
dynamics- are described by a measure in the So-matrix 
space which, through the assumption of ergodicity, gives 
the probability of finding So in a given volume element as 
the energy E changes and So wanders through that space 
4]. In the absence of direct processes, or short paths, So 
is taken to have a uniform distribution, given by the in- 
variant measure for the symmetry class in question. The 
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invariant measure, denoted by dfip(So), is denned by its 
invariance under the s ym metry operation relevant to that 
universality class |2f|[2]|; i.e., 



d/J.p(S ) = dfi (U o S o V o ), 



(2.3) 



where Uq and Vq are arbitrary but fixed unitary matrices 
in the unitary case, while Vq = Uq in the orthogonal case. 
This defines the circular unitary (orthogonal) ensemble, 
CUE (COE) for /3 = 2 (fl = 1). 
The distribution 



dPp{S)= P p{S)d^{S) 



(2.4) 



of the resulting scattering matrix S of Eq. 1(2. 2|1 for the 
full system consisting of the cavity plus the absorbing 
mirror can in principle be calculated for a given S a , as- 
suming that So is distributed according to the invariant 
measure, Eq. ((2.3(1 . This procedure is particularly well 
suited for numerical simulations. Indeed, it is used in 
Sec. II VI to study numerically the case of one open chan- 
nel (N = 1) and one absorbing channel (N a = 1). 

One property of pp(S) is immediately obvious: the av- 
erage of S vanishes, i.e., (S) = 0, as we can see by ex- 
panding the right-hand side of Eq. 1)2.2(1 and averaging 
term by term. 

Although the explicit expression for pp{S) is not known 
for an arbitrary number of channels, the case of one open 
channel [N — 1) and one absorbing channel (N a = 1) can 
be solved analytically: this is done in Sec. IIVI exploit- 
ing the relation with the parasitic-channel model that is 
explained in the next section. 



III. RELATION BETWEEN THE 
SURFACE-ABSORPTION MODEL AND THE 
PARASITIC-CHANNEL MODEL 

In the parasitic-channel model described in the Intro- 
duction and sketched in Fig. [2] the cavity is attached 
to a waveguide supporting N physical channels, while 
absorption is simulated by attaching an additional, fic- 
titious, waveguide supporting N p equivalent "parasitic" 
channels, each one with a coupling given by a transmis- 
sion coefficient T p @, [13 . 

We denote by S the scattering matrix that describes 
the whole system shown in Fig. [21 including the fictitious 
waveguide and the (non-absorbing) barrier inside it; S is 
of dimension N + N v and has the structure 



5 = 



S Sip 
Sp i Spp 



(3.1) 



where the indices 1 and p denote the set of N physical and 
N p fictitious channels, respectively; it is unitary for (3 — 2 
and unitary-symmetric for /3 = 1 . The N x N submatrix 
S describes the cavity in the presence of absorption, as 
one would observe it in an experiment in which only the 




FIG. 2: In the parasitic-channel model, a ballistic chaotic 
cavity connected to a waveguide with N propagating modes 
is also connected to an additional, fictitious, waveguide sup- 
porting N p "parasitic" channels. Inside the latter waveguide 
there is a (non-absorbing) barrier described by the scattering 
matrix S&. Only the N x N submatrix S of the total scat- 
tering matrix of this system is accessible to experiment; it is 
subunitary and is taken as a model for a real cavity in the 
presence of absorption. 



physical waveguide were present; S is subunitary and will 
be referred to as the "physical" S matrix. 

The tunnel barrier in the fictitious waveguide in Fig. 
13 described by the 2N P x 2N p unitary scattering matrix 
S b 



(3.2) 



implies an imperfect coupling of the parasitic channels 
to the cavity; it can also be interpreted as giving rise to 
direct processes in a system described by S [IJ. Then, 
for a cavity whose classical dynamics is chaotic, S is dis- 
tributed according to the Poisson kernel |28l 129. 30] 

(/3N T +2-0)/2 



det 



(t NT -(S)(S) 



det (1^ - 5(5) t 



f3N T +2-P 



-dnp{S), 



(3.3) 

where Nt = N + N p , dfj,p{S) is the invariant measure for 
5 defined as in Eq. 12. 3D . assumed to be normalized to 
unity, and 



(5) 



N _ 

o 







(3.4) 



The probability distribution of the physical scattering 
matrix 5 that we are interested in has to be calculated 
from Eq. I|3.3|) . This is how Brouwer and Beenakker 
showed that, in the limit N p — > oo and T p — > 0, with 
the product 7 p = N p T p being kept fixed and interpreted 
as the absorption strength, this model describes volume 
absorption, in the sense that the problem represented by 
5 is equivalent to one in which all the energy levels of the 
closed cavity acquire a fixed imaginary part, associated 
precisely with that absorption strength 0, |Tl| . 

Let us now think of the full system shown in Fig. [21 
as the combination of the cavity connected to the two 
waveguides (in the absence of the barrier) and described 
by the (N + N p ) x (N + N p ) scattering matrix 



So 



(3.5) 
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plus the barrier described by Sb- We shall not take, in 
what follows, the limit mentioned in the previous para- 
graph, i.e., Np will be kept finite. As we shall see, it 
will also be convenient, for the description of the barrier, 
to drop the "equivalent-channel" assumption normally 
made in this context [Eq. l|3.2|l ] and consider a barrier 
matrix Sb with the general structure 



S h = 



n t' b 



The total scattering matrix S is then given by 

1 



S = r B + t' B 



Iw+Wp — S r' B 



Sots 



(3.6) 



(3.7) 



Here we have introduced the (N+N p ) x (N+N p ) matrices 



t B 



On 

n 

Ijv 

t b 








(3.8a) 
(3.8b) 



which denote the barrier reflection and transmission ma- 
trices associated with the two waveguides in Fig. [2] [no- 
tice that in the present model the left waveguide has no 
barrier; hence the matrices Oat and 1at in Eq. Ij3.8|l : in 
a more complete model there might be a barrier in the 
left waveguide too (see comment at the end of Sec. II V|) ]. 
The physical S matrix S of Eq. I|3.1fl can then be found 
from Eq. I|3.7|l to be 

1 



S = f o + t' Q 



N„ 



r 'b f 



r' h t 



b L o 



(3.9) 



Comparing Eqs. (|2.2(l and l|3.9(l we see that the model 
for surface absorption of Sec. [H]and the parasitic-channel 
model of the present section can be made equivalent, in 
the sense that S — 5, choosing the number of parasitic 
channels in the latter to coincide with the number of ab- 
sorbing channels in the former, i.e., N p — N a , choosing 
the two cavities to be identical, i.e., So = So, and the bar- 
rier reflection matrix r' b of the present section to coincide 
with the absorption matrix of Sec. [n] i.e., r' b = S a - 

Within the equivalent-channel assumption, equiva- 
lence of the two models would require S a — yj\ — T p t m p ■ 
If in our surface-absorption model we were to take the 
limit that was described right after Eq. I|3.4[) in connec- 
tion with the parasitic-channel model, we would obtain 
the same answer as for volume absorption by identifying 
7 P with 7 a = trP Q where P a = t Na - S\S a . 

The statistical distribution of the scattering matrix S 
of Eq. H3.7JI is of course given by Poisson's kernel of Eq. 
B3.3J) . where the average S matrix (S) is now given by 
(S) = r B [see Eq. (l3~%a|) ]. instead of Eq. (|3~4"|) . Using 
the identification described right after Eq. (|3.9|l , the cal- 
culation of the resulting distribution of S also gives the 
distribution of S, Eq. (|2.2J) . for the model of the previous 
section. This procedure will be taken advantage of in the 
next section, in order to analyze the case N = 1, N a = 1. 



IV. ANALYTICAL SOLUTION OF THE 
SURFACE-ABSORPTION MODEL FOR THE 
CASE N = 1, N a = 1 

In this section we calculate, within the surface- 
absorption model of Sec. [H] the distribution of the scat- 
tering matrix S of Eq. (|2.2|l which describes a chaotic 
cavity perfectly connected to a waveguide supporting one 
channel (N = 1), and in the presence of one absorption 
channel (N a = 1); S is thus a complex number that we 
parametrize as 



S = VRe l 



(4.1) 



The absorbing mirror is also described by a subunitary 
matrix S a = \ZR a e lSa . 

We can take advantage of the equivalence with the 
parasitic-channel model discussed in the previous sec- 
tion, with the identification: N p = N a = 1 and T p — 

P a = 1 — R a , a = 0. For the present case, S of Eq. I|3.1() 
is a 2 x 2 unitary matrix for (3 — 2, with the additional 
condition of symmetry for (3=1. It can be parametrized 
in a polar representation as [3^, |3^| 



S = 



vr 



're 



re 



+4>') 



v 7 ! 



(4.2) 



where < <f>, ip, (/>', tp' < 2ir and < r < 1; for j3 = 1 we 
have the restrictions <p' — cf> and tp 1 — tp. The 11 matrix 
element of S is the physical S = S = -VT _7 Te i ^+^'), 
so that the parameters of S appearing in (|4.1|l are given 
by 



R = 1- 



The matrix S is distributed according to Eq. 
N T = N + N p 



2 and 

(S) = 





IRa 



(4.3a) 
(4.3b) 

with 



(4.4) 



In the polar representation of Eq. 14. 2|) . the invariant 
measure dnp(S) needed in Eq. I|3.3|l is given by [31111 

- dr dtpdip f d(j)' dip'Y^ 1 , AK ^ 



The quantity 
(|3~3|) . i.e., 



(2 v /t) 2 - /3 2tt 2tt V 2tt 2tt 



det (l N+Np -S(S)l 



appearing in Eq. 



det 



(t 2 -s(sy 



= i + (i - r)R a - 2VT^VRa cos(v> + 

(4.6) 

is independent of the angles <f>, <j)' and hence of the phase 9 
of S, which is thus uniformly distributed. The normalized 
probability distribution of S can thus be written as 

d6 



dPp(S) = P P {R) dR 



2tt' 



(4.7) 



What remains to be calculated is thus the probability 
distribution of the reflection coefficient R, Pp(R). Refs. 
[8?il | and Q give the distribution of r for certain particular 
cases of the average S matrix, of which Eq. 1)4. 4f> is a 
particular case. In Eqs. (5.8) and (5.11) of Ref. [3(j 
(or Eqs. (6.81) and (6.85) of Ref. 0), we make the 
substitution 1 — T = R, together with X = 0, Y = \/R^ 
in (5.8) and X = y/R^, Y = in (5.11) of Ref. [13, (or 
Eqs. (6.81) and (6.85) of Ref. Q, respectively), with the 
result 



^2 OR) 



(l-R a ) 2 {l + R a R) 

(1 - RaR) 3 



for (3 = 2, and 



(1 - Raf /2 

2V 7 ! - !? 



2^(3/2,3/2;^) 



(4.8) 



(4.9) 



for fi = 1, 2-Fi being a hypergeometric function [361 ]. 

In principle, the distribution of S could also be ob- 
tained from its definition, Eq. 1)2. 2JI . and the fact that So 
is distributed according to the invariant measure. This 
procedure is especially suited for numerical simulations. 
The absorbing mirror is described by the lxl subuni- 
tary scattering matrix S a = V R a e l0a , where the phase 
9 a and the reflection coefficient R a are fixed; < R a < 1 
and for simplicity we may take, as above, 9 a = 0. The 
lack of unitarity of S a is given by P a = 1 — R a . The 
matrix So of Eq. 1)2.1)1 is a 2 x 2 unitary matrix that 
can be parametrized in a polar representation as in Eq. 
1)4.2(1 . with the various parameters now having an index 
and the range of variation: < 4>' , tp' , , ip' ' < 2tt and 
< To < 1; for (3 = 1 we have the restrictions cf>' Q = <f>o 
and ip' — ipQ. For a chaotic cavity, So is distributed ac- 
cording to the invariant measure which in terms of the 
polar parameters is given by 



dr dcj) dipp / d<t>' Q dip' Q 
(2^) 2 "' 3 2tt 2vr \ 2tt 2tt 



8-1 



r^zp . (4.10) 



The distribution of S can then be obtained from Eq. 
(I2~2l . 

In Fig. |21 we compare the analytical results of Eqs. 
(|4.8f) and 1)4.9)1 with the results of a numerical simulation 
obtained from the ensemble of Eq. 1)4.10)1 and the relation 
(12.21) between S and So: the agreement is excellent. 

We should point out that in real experiments the 
waveguide is not perfectly coupled to the cavity, as has 
been shown in Ref. [l5| . thus giving rise to direct pro- 
cesses. In such cases, the distribution of S must be mod- 
ified by the procedure of Ref. [2lJ • 



V. SUMMARY AND CONCLUSIONS 

We have presented a model to describe the statistical 
properties of wave scattering in a ballistic chaotic cav- 
ity in the presence of surface absorption: it is suggested 
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FIG. 3: Distribution of R of Eq. J£TJ for a) /3 = 2 and 
b) fi = 1 symmetries, for one open channel, N = 1, one 
absorbing channel, N a = 1, and R a = 0.1, 0.5, and 0.9. The 
continuous lines are the analytical results given by Eqs. lETgt 
for (3 = 2 and 14.911 for fi = 1. The numerical simulation was 
performed as explained in the text using 10 4 samples; the 
different values of R a axe distinguished by the symbols also 
shown in the figure. The agreement between theoretical and 
numerical results is excellent. 



that this type of absorption be realized in the labora- 
tory by attaching one or several "absorbing patches" to 
the inner wall of the cavity. The model simulates one 
such absorbing patch by means of an absorbing mirror, 
whose physical characteristics, i.e., its cross section and 
its absorption properties, are, in principle, experimen- 
tally measurable. 

We compare our model with the parasitic-channel 
model introduced by other authors and find the condi- 
tions under which the two are equivalent. It is found 
in the literature that in the limit of a large number 
of parasitic channels weakly coupled to the cavity, the 
parasitic-channel model describes absorption uniformly 
distributed over the volume of the cavity; in this limit, 
volume and surface absorption give the same results. 

Finally, we show that our model is analytically solvable 
for the case of one open channel, N = 1, and one absorb- 
ing channel, N a = 1, for arbitrary absorption strength. 
The results of numerical simulations for this case are in 
excellent agreement with theory. In principle, this prob- 
lem is amenable to experimental observation. 
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